Quadruple inversion recovery for quantitative contrast-enhanced black blood imaging

ABSTRACT

A contrast enhancement (CE) agent is infused into blood flowing through a site that is to be imaged with magnetic resonance imaging (MRI). Two double inversion procedures are carried out, forming a quadruple inversion recovery (QIR) pulse sequence. Each double inversion procedure comprises a non-selective and slice-selective inversion RF pulse. The first double inversion procedure is followed by a first predefined inversion delay period, TI 1 , and the second procedure by a second predefined inversion delay period, TI 2 . A black-blood image can thus be produced in which blood appears consistently black and tissues surrounding the blood, such as a vessel wall, heart, atherosclerotic plaque, or thrombus, are clearly visible. Unlike the prior art black-blood imaging technique, the QIR method does not require a precise knowledge of the T 1  of the blood carrying the CE agent in order to suppress the signal and artifacts caused by the blood flowing through the site.

RELATED APPLICATIONS

[0001] This application is based on prior copending provisional patent application, Serial No. 60/435,835, filed on Dec. 19, 2002, the benefit of the filing date of which is hereby claimed under 35 U.S.C. §119(e).

GOVERNMENT RIGHTS

[0002] This invention was made with government support under Grant No. RO1-HL56874, awarded by NIH-NIHBI, and the federal government may have certain rights in this invention.

FIELD OF THE INVENTION

[0003] This invention generally pertains to magnetic resonance imaging (MRI) and more specifically, pertans to a method, a pulse sequence, and a system that make use of a plurality of radiofrequency (RF) inversion pulses grouped into double-inversion procedures in order to suppress the signal from flowing blood having a variety of T₁ relaxation times caused by the presence of a contrast-enhancing agent.

BACKGROUND OF THE INVENTION

[0004] Black-blood imaging is a technique developed to improve the visualization of blood vessels and the heart in magnetic resonance images so as to more clearly reveal pathologic tissues and morphologic abnormalities, which appear close to the blood-wall interface. (See Edelman R. R., Chien D., and Kim D., “Fast Selective Black Blood MR Imaging,” Radiology 1991; 181:655-660.) These improvements in visualization include the elimination of image artifacts arising from flowing blood and a reduction in the residual blood signal, which may mask pathologic conditions such as atherosclerotic plaque or thrombus formation. The combination of black-blood imaging with contrast enhancement (CE) offers a high potential for various cardiovascular applications, and in particular, for improved high-resolution MRI of atherosclerotic plaque. (See Yuan C. et al., “Contrast-Enhanced High Resolution MRI for Atherosclerotic Carotid Artery Tissue Characterization,” Journal of Magnetic Resonance Imaging, 2002; 15:62-67.) The ability to provide high resolution images of plaque in the carotid arteries and other vessels is becoming increasingly important in assessing the potential vulnerability and risk of patients to stroke and other cardiovascular diseases. However, due to the significant shortening of T₁ in blood and the variability of T₁ caused by the use of a contrast enhancement (CE) agent, a traditional method for achieving black-blood imaging, such as double inversion-recovery (DIR) cannot guarantee effective blood suppression on post-contrast images. This issue is especially important for quantitative image analysis, such as the calculation of contrast enhancement and morphological measurements. The crucial problem of DIR is that it is necessary to know the relaxation time, T₁, of blood in order to calculate the inversion time (TI). The DIR method can effectively suppress the blood signal if and only if TI is determined so as to allow the magnetization of blood to approach zero. However, the use of CE agent decreases the T₁ parameter of blood by an unpredictable amount. Uncertainties in the effect of CE on T₁ arise from various factors related to the injection technique, flow dynamics, concentration, and time between injection and imaging. Determining the proper TI is therefore uncertain. Accordingly, it would be critical to employ an MRI technique that provides the benefits of DIR in achieving black-blood images, but is much less sensitive to variations in the T₁ parameter of blood due to application of a CE agent.

SUMMARY OF THE INVENTION

[0005] As explained above, the current approach used for black-blood imaging in MRI of the heart and blood vessels has typically applied the DIR method. However, DIR only provides efficient suppression of the blood signal within a narrow range of T₁ values at any selected TI. The present invention addresses this issue by suppressing the signal and artifacts from flowing blood in MRI of the heart and blood vessels over a wide range of T₁ values, so that the changes in T₁ caused by the infusion of a CE agent are not problematic in achieving high-quality black blood images.

[0006] More specifically, one aspect of the present invention is directed to a method for suppression of a signal from flowing blood in MRI, so that tissue surrounding the blood is clearly visible in the image, while the blood is consistently black. The method comprises the step of applying two (or more) double-inversion procedures before signal acquisition. The method starts when a first double-inversion procedure is applied. After waiting for a first predefined inversion delay period, a second double-inversion procedure is carried out. Each double-inversion procedure comprises a non-selective RF pulse applied to provide a non-selective inversion of the magnetization in the entire object and a slice-selective pulse that is targeted to invert the magnetization in a slice, a group of slices, or a three-dimensional (3-D) slab. After waiting for a second predefined inversion delay period, a signal is acquired for processing to produce an image in which the blood flowing through the slice is black, in contrast to tissue comprising the blood vessel wall, heart, and any pathologic tissue such as atherosclerotic plaque or thrombus that may have been formed inside the blood vessel. The first and second predefined inversion delay periods are selectively determined so that effects of the flowing blood are suppressed in the image over a wide range of T₁ values, which may occur in the blood both in the presence and in the absence of the CE agent. This result is achieved when a longitudinal magnetization of the blood flowing through the site approaches zero at a time when the signal is acquired for imaging.

[0007] The longitudinal magnetization of the blood varies as a combination of exponential functions of several parameters. These parameters include the first predefined inversion delay period, the second predefined inversion delay period, a repetition time, TR, and a current T₁ value of the blood, which depends on the concentration of the CE agent infused. The CE agent reduces the T₁ value for blood in an unpredictable manner. However, the residual magnetization of the blood is negligibly small for a wide range of T₁ after application of the sequence of RF pulses in accord with the present invention. This range covers all possible T₁ variability. Thus, the method produces consistently black-blood images independently of the presence and concentration of the CE agent in the blood.

[0008] In another aspect of the invention, a method for black-blood contrast enhanced MRI is defined. The method starts with the step of infusing the CE agent into the blood. The CE agent used preferably includes gadolinium derivatives, although other suitable CE agents may instead be used. After waiting a sufficient time, blood carrying a portion of the CE agent flows through the site. At that time, the images are acquired using the method of the present invention, which includes two double-inversion procedures, each of which is followed by one of two predefined time periods. As was noted above, this method effectively suppresses the blood signal. Thus, tissues surrounding blood are clearly visible after CE without any problems associated with flow artifacts and a residual blood signal.

[0009] In addition, the black-blood image can be obtained at the site before infusion of the CE agent. Essentially, this image guarantee consistent suppression of the signal from blood using the same parameters that are used in the present invention for post-CE imaging. Because images before and after CE have been obtained using the identical method, any changes of signal intensity are caused only by the penetration of the CE agent into tissue. This information disclosed by the resulting images is used for the characterization of pathologic tissues in the cardiovascular system, such as atherosclerotic plaques, thrombi, myocardial infarction, and other conditions.

[0010] Another aspect of the invention is directed to a pulse sequence for quadruple inversion-recovery imaging. The pulse sequence includes two double inversion RF pulse pairs, each of which is followed by one of two inversion delays. Each double inversion pulse pair includes two inversion RF pulses, one of which is non-selective and the other of which is slice-selective or slab-selective. The pulse sequence also includes all RF and magnetic field gradient pulses needed to generate, encode, and acquire an MR signal. Optionally, the sequence may contain more than two double inversion blocks in order to improve blood suppression effect.

[0011] Yet another aspect of the invention is directed to a system for quadruple inversion-recovery data acquisition. The system includes a computer having a processor and a memory in which a set of machine instructions are stored. The machine instructions are executed by the processor. The computer is connected to an MRI apparatus and execution of the machine instructions causes the MRI apparatus to generate the quadruple inversion recovery (QIR) pulse sequence. The MRI apparatus is thus controlled by the computer to produce images using magnetic field and RF pulses as discussed above. The machine instructions also determine appropriate inversion delays to achieve black-blood imaging m which the blood signal is substantially zero. The execution of the QIR pulse sequence can be performed with a TR or synchronized with the cardiac rhythm of the living object (human or animal).

BRIEF DESCRIPTION OF THE DRAWING FIGURES

[0012] The foregoing aspects and many of the attendant advantages of this invention will become more readily appreciated as the same becomes better understood by reference to the following detailed description, when taken in conjunction with the accompanying drawings, wherein:

[0013]FIG. 1 is a graph for a simulation showing the magnetization of blood as a function of T₁ for DIR and QIR sequences, at a TR equal to 800 ms;

[0014]FIG. 2 is a graph of normalized signal intensity in relation to pre-CE and post-CE and TI for a DIR-fast-spin echo (FSE) sequence (with IR equal to 800 ms) measured using fresh bovine muscle;

[0015]FIG. 3 illustrates T₁-weighted QIR-FSE MRI images, pre-CE (top) and post-CE (bottom), of a carotid atherosclerosis patient;

[0016]FIG. 4 is a timing diagram of the FSE-QIR sequence; and

[0017]FIG. 5 is a block diagram showing a computer coupled to an MRI apparatus, for practicing the present invention.

DESCRIPTION OF THE PREFERRED EMBODIMENT

[0018] Method

[0019] Sequence: In the present invention, suppression of the flowing blood signal and artifacts is achieved using a quadruple inversion recovery (QIR) pulse sequence. The QIR pulse sequence includes two double inversion pulse pairs, the first double inversion pulse pair being followed by an inversion delay period, TI₁, and the second double inversion pulse pair being followed by an inversion delay period, TI₂. The double inversion procedure is known from its use in a double inversion-recovery (DIR) pulse sequence. Use of the DIR sequence for achieving black-blood imaging is generally well known to those of ordinary skill in the art of MRI. It is also known that the DIR method provides efficient suppression of the blood signal in a narrow range of T₁ values, which makes the use of DIR problematic in the presence of contrast enhancing agents, since such agents strongly affect the T₁ of blood. However, in accord with the present invention, by performing a plurality of successive double inversion procedures, each followed by an appropriate inversion delay period, black-blood imaging can be achieved without requiring a precise knowledge of T₁ in blood after infusion of the contrast enhancement agent. Although initially tested as a QIR sequence including two double inversions, it is contemplated that the present invention may be implemented as a sequence of more than two double inversions, each followed by an appropriate inversion delay period.

[0020]FIG. 4 illustrates the steps employed in the QIR pulse sequence. As shown therein, a first double inversion 10 is carried out, followed by a first inversion delay period, TI₁, which is identified by reference number 12. A second double inversion 14 is then implemented, followed by a second inversion delay period, TI₂, which is identified by a reference number 16. A signal produced by the tissue being imaged is then acquired, for example, using an FSE pulse sequence, as indicated in a block 18. Generally, any two-dimensional (2-D) or 3-D acquisition sequence can be used to produce signal excitation and spatial encoding instead of the FSE pulse sequence.

[0021] In each double inversion procedure, a slice-selective inversion radio frequency (RF) pulse is applied to the site to be imaged, where the blood flow occurs (e.g., blood vessels or the heart) immediately after a nonselective inversion RF pulse has been applied. The inversion pulses also can be applied in the reverse time order with the same effect. The nonselective pulse inverts the magnetization of the entire object placed inside a transmit coil of the scanner, which is generally the part of the patient's body covered by the coil and includes a substantial volume of circulating blood. The slice selective pulse re-inverts the magnetization of the imaged slice to its initial direction. Also, it should be noted that adiabatic pulses can be used for either the selective or the nonselective inversion in order to reduce phase errors of magnetization across the object. The thickness of the re-inverted slice is generally greater (1.5-3 times) than the thickness of the imaged slice in order to avoid displacement effects caused by breathing and heart pulsation. The slice selective pulse can be applied to a single slice, or a group of slices, or a three-dimensional slab depending on the acquisition method (2-D single- or multi-slice, or 3-D). If the first and second inversion delay periods, TI₁ and TI₂ are properly selected, the longitudinal blood magnetization, M_(z), should approach zero after the second inversion delay period, TI₂. While the first and second inversion delay periods can be mathematically calculated, it is also possible to empirically determine each of the inversion delay periods using graphical methods.

[0022] In an initial test of the present invention that confirmed its ability to provide improved black-blood imaging after an injection of the contrast enhancement agent, the QIR pulse sequence was implemented on a 1.5T MR scanner (GE Signa EchoSpeed™) in combination with FSE readout sequence. It will be understood that other suitable image acquisition techniques can be implemented to obtain the signal from a site produced using MRI. The QIR pulse sequence can be implemented with two timing options. In the first option, the repetition time (TR) or the sequence is fixed and defined by the operator. In the second option, the execution of the sequence is synchronized with the cardiac rhythm of the patient so that the repetition time is equal to one or more R-to-R intervals. This cardiac-gated acquisition mode is necessary for imaging of anatomic locations affected by heart pulsation such as the heart itself and the aorta. In the tests demonstrating the efficacy of the present invention, the QIR pulse sequence was implemented with adiabatic hyperbolic-secant RF pulses applied for slice-selective inversion and rectangular RF pulses applied for non-selective inversion. It further will be understood that other pulse shapes can be used in the QIR pulse sequence to achieve spatially selective and nonselective inversion of magnetization.

[0023] Theory: Longitudinal magnetization of blood outside an imaging slice under the action of the QIR pulse sequence is described by the following equation: $M_{z} = \frac{1 - {\exp \left( {- \frac{TR}{T_{1}}} \right)} - {2\quad {\exp \left( \frac{- {TI}_{2}}{T_{1}} \right)}} + {2\quad {\exp \left( \frac{- \left( {{TI}_{2} + {TI}_{1}} \right)}{T_{1}} \right)}}}{1 - {\exp \left( \frac{- {TR}}{T_{1}} \right)}}$

[0024] which has two zero solutions at predefined TR, TI₁, and TI₂ (see the example shown in FIG. 1). In addition, longitudinal magnetization of blood outside an imaging slice is negligibly small over a wide range of T₁, as seen in FIG. 1, which accounts for the feasibility of suppressing blood signal with a variety of T₁ values without knowledge of the actual T₁.

[0025] MRI Experiments: The effect of DIR and QIR pulse sequences on a steady tissue signal was studied ex vivo using a sample of fresh bovine muscle. Subsequently, to test the QIR sequence in vivo, high-resolution T₁-weighted images of carotid arteries were obtained from two carotid atherosclerosis patients (before and after contrast enhancement agent administration) and two healthy volunteers (who did not receive any contrast enhancement injection). Informed consent was obtained from all participants. The protocol included DIR (inversion time, TI=325 ms) and QIR (TI₁/TI₂=375/125 ms) scans, with the following parameters of the FSE acquisition sequence: TR=800 ms, TE=11 ms, echo train length=8, two signal averages, matrix 256×256, field of view 13×7 cm, and a slice thickness of 2 mm. Post-contrast scans were acquired 5-15 min. after an injection of 20 ml of a gadolinium-based contrast agent (Omniscan™).

[0026] Results of the Experiments

[0027] As was shown by simulations and the ex vivo experiments, the following limitations of the DIR method hampers quantitative contrast enhanced studies:

[0028] 1. The DIR method provides efficient suppression of the blood signal within only a narrow range of T₁ values at any chosen TI (see FIG. 1, which illustrates two plots of blood magnetization for the DIR sequence, with TI=325 ms, and with TI=175 ms). Thus, DIR is prone to the blood flow signal and artifacts that are caused by variations of T₁ in blood after a contrast enhancement agent administration.

[0029] 2. DIR is characterized by strong (up to 20% of relative change) dependence of a steady tissue signal on TI (as shown in FIG. 2). It has been demonstrated that this dependence is determined mainly by the magnetization transfer (MT) effect induced by the pair of inversion pulses used in the DIR procedure. While the single DIR procedure should be applied with different TI for pre- and post-contrast enhancement infusion scans to obtain satisfactory blood signal suppression, it may cause a misinterpretation of the contrast enhancement in weakly enhanced tissues.

[0030] Features of the QIR sequence. The QIR sequence with properly chosen TI₁ and TI₂ provides efficient blood signal suppression over a relatively wide range of TI values (as shown by the solid line for the QIR sequence in FIG. 1). In the present protocol, the parameters of the QIR-FSE sequence were adjusted to guarantee at least 90% blood signal suppression, starting from an expected T₁=200 ms. The predefined values for TI₁ and TI₂ used were sufficient for obtaining perfect lumen visualization in both arterial and venous vessels (as shown in FIG. 3). FIG. 3 includes images made with the QIR sequence in accord with the present invention for both pre-contrast enhancement agent injection (top image) and post-contrast enhancement agent injection (bottom image).

[0031] Note that unlike DIR, QIR does not require a precise knowledge of T₁ for blood after a contrast enhancement agent injection, since the corresponding suppression curve is quite flat (see the solid line curve for the QIR pulse sequence shown in FIG. 1). While the QIR sequence is subjected to a decrease of steady tissue signal due to MT saturation, similar to that of DIR, both pre-, and post-contrast QIR images contain the same amount of MT effect, because they were obtained with the identical QIR sequence. As a result, it is possible to attribute quantitative changes of signal intensity to the influence of a contrast enhancement agent only. Although the saturation of steady tissues by the QIR sequence appeared to be slightly higher than that for DIR preparation (3-9% of further signal decrease), this effect was too small to affect the overall image quality.

[0032] System for Implementing the Present Invention

[0033]FIG. 5 schematically illustrates a system for implementing the present invention that includes a generally conventional MRI apparatus 30 that is controlled by a computer 32. Computer 32 may be a generally conventional personal computer (PC) or a dedicated controller specifically intended for controlling MRI apparatus 30. Although not specifically shown, MRI apparatus 30 includes a magnet to create permanent magnetic field, a plurality of gradient coils to produce spatial variations of magnetic field, and an RF transceiver and receiver systems to transmit and receive RF signals to and from a plurality of RF coils, as will be well known to those of ordinary skill in the art of MRI. Accordingly, details of the MRI apparatus need not be specifically illustrated or discussed herein.

[0034] Computer 32 is coupled to a display 34 which is used for displaying MRI images to an operator. Included within computer 32 is a processor 36. A memory 38 (with both read only memory (ROM) and random access memory (RAM), a storage 40 (such as a hard drive or other non-volatile data storage device) for storage of data and software programs, an interface 44, and a compact disk (CD) drive 46 are coupled to processor 36 through a bus 42. CD drive 46 can play a CD 48 on which machine instructions are stored for implementing the present invention and other software modules and programs that may be run by computer 32. The machine instructions are loaded into memory 38 before being executed by processor 36 to carry out the steps of the present invention.

[0035] Operation of MRI apparatus 30 is controlled by computer 32 when processor 36 executes the machine instructions stored in memory 38. These machine instructions cause the processor to implement the sequence using two or more double inversion procedures, each followed by an inversion delay period, as discussed above. The resulting image is displayed on display 34 using data stored in storage 40 to enable further processing. It is also noted that the machine instructions will cause processor 36 to determine the appropriate inversion delay periods following each double inversion procedure so as to minimize a deviation of the blood magnetization from zero over a relatively wide range of T₁.

[0036] Conclusion

[0037] The QIR pulse sequence in accord with the present invention allows high-quality black-blood images to be obtained both before and after administration of the gadolinium-based contrast enhancement agent. The method of the present invention thus provides reliable quantitative interpretation of a contrast enhancement in tissues and can be employed for MRI in a wide spectrum of cardiovascular applications.

[0038] Although the present invention has been described in connection with the preferred form of practicing it, those of ordinary skill in the art will understand that many modifications can be made thereto within the scope of the claims that follow. Accordingly, it is not intended that the scope of the invention in any way be limited by the above description, but instead be determined entirely by reference to the claims that follow. 

The invention in which an exclusive right is claimed is defined by the following:
 1. A method for suppressing a signal from flowing blood in a signal used for magnetic resonance imaging (MRI) of a site in a subject, comprising the steps of: (a) generating a sequence of radiofrequency (RF) pulses, a total number of the said RF pulses being an even number at least equal to four; (b) applying said RF pulses to the subject so that the said RF pulses are grouped into a plurality of double inversion procedures, each double inversion procedure including a non-selective RF pulse targeted to the subject and a slice-selective RF pulse that is targeted to one of a slice, a group of slices, and a three-dimensional slab at the site, wherein one of the non-selective RF pulse and the slice-selective RF pulse immediately follows the other of the non-selective RF pulse and the slice-selective RF pulse; (c) waiting a predefined time period after each of the said double inversion procedures; and (d) after the predefined time period following a last of the said double inversion procedures, acquiring a signal that is processed to produce an image of the site, said signal from the flowing blood being substantially suppressed in the signal used to produce the image of the site.
 2. The method of claim 1, wherein the predefined time periods are determined so that the longitudinal magnetization of flowing blood is substantially zero when the signal used to produce the image is acquired.
 3. The method of claim 1, further comprising the step of infusing a contrast-enhancing agent in the subject before generating the sequence of RF pulses, wherein the predefined time periods are determined so that the signal from flowing blood is sufficiently suppressed at any value of a T₁ parameter affected by the infusion of the contrast-enhancing agent.
 4. The method of claim 1, wherein the sequence of RF pulses comprise a first double inversion procedure and a second double inversion procedure, the first double inversion procedure being followed by a first predefined time period, TI₁, and a second double inversion procedure being followed by a second predefined time period, TI₂.
 5. The method of claim 4, wherein the longitudinal magnetization of the blood is defined as a function of: (a) a repetition time, TR; (b) a T₁ parameter for the flowing blood; (c) the first predefined time period TI₁; and (d) the second predefined time period TI₂.
 6. The method of claim 5, wherein the first and the second predefined time periods, TI₁ and TI₂, are determined so that the longitudinal magnetization of the flowing blood is substantially zero over a wide range of the T₁ parameter.
 7. A method for suppressing the signal from flowing fluid in magnetic resonance imaging (MRI) after infusion of a contrast agent, comprising the steps of: (a) carrying out a first double inversion procedure that includes applying in rapid succession and in either order a non-selective inversion RF pulse and a slice-selective inversion RF pulse, said slice-selective inversion RF pulse being targeted at one of a slice, a group of slices, and a three-dimensional slab; (b) waiting a first predefined inversion delay period; (c) carrying out a second double inversion procedure that includes again applying in rapid succession and in either order a non-selective inversion RF pulse and a slice-selective inversion RF pulse targeted at said one of the slice, the group of slices, and the three-dimensional slab; (d) waiting a second predefined inversion delay period; and (e) acquiring a signal for processing to produce the image in which the signal from flowing fluid appears suppressed and a tissue surrounding the fluid is consistently visible.
 8. The method of claim 7, wherein the fluid comprises blood, and the tissue surrounding the fluid comprises one of a wall of a blood vessel and a heart.
 9. The method of claim 7, wherein the first predefined inversion delay period and the second predefined inversion delay period are selected so that a longitudinal magnetization of the fluid is substantially zero when the signal is acquired for processing to produce the image.
 10. The method of claim 7, wherein a longitudinal magnetization of the fluid is defined as a function of: (a) a repetition time, TR; (b) a T₁ parameter for the fluid; (c) the first predefined inversion delay period; and (d) the second predefined inversion delay period.
 11. The method of claim 7, wherein the contrast enhancement agent includes but not limited to a gadolinium derivative.
 12. The method of claim 7, wherein the use of the first and the second double inversion procedures suppresses a contribution of the fluid to the image over a wide range of T₁, for the fluid with the contrast enhancement agent infused therein.
 13. A memory medium storing machine instructions for carrying out the steps of claim
 1. 14. A method for contrast-enhanced black-blood magnetic resonance imaging (MRI) so that tissue surrounding blood is clearly visible and the blood is consistently black in an image of a site, comprising the steps of: (a) infusing a contrast enhancement agent into the blood so that the blood carrying a portion of the contrast enhancement agent is flowing through the site to be imaged; (b) carrying out a first double inversion procedure; (c) waiting a first predefined inversion delay period; (d) carrying out a second double inversion procedure, each of the first and the second double inversion procedures including the steps of applying in rapid succession an in either order a non-selective inversion RF pulse and a slice-selective inversion RF pulse, said slice-selective inversion RF pulse being targeted at one of a slice, a group of slices, and a three-dimensional slab at the site; (e) waiting a second predefined inversion delay period; and (f) acquiring a signal for processing to produce an image, in which the flowing blood is consistently black.
 15. The method of claim 14, wherein the first and second predefined inversion delay periods are selectively determined so that a signal from the flowing blood is suppressed in the image over a relatively wide range of T₁ values for the blood carrying the portion of the contrast enhancement agent.
 16. The method of claim 14, wherein a longitudinal magnetization of the flowing blood is substantially zero at a time when the signal is acquired.
 17. The method of claim 14, wherein said longitudinal magnetization varies as a functions of the first predefined inversion delay period, the second predefined inversion delay period, a repetition time (TR), and a T₁ value of the blood carrying the portion of the contrast enhancement agent.
 18. The method of claim 14, wherein the contrast enhancement agent comprises a gadolinium derivative.
 19. The method of claim 14, further comprising the step of applying the first double inversion and the second double inversion, each followed by the predefined inversion delay periods, so that a signal is acquired for producing a non-enhanced image obtained before the contrast enhancement agent is infused into the blood.
 20. The method of claim 19, wherein the said non-enhanced image is obtained with inversion delay periods and a repetition time (TR) that are the same as the inversion delay periods and the repetition time (TR) used when producing the image after the contrast enhancement agent was infused.
 21. The method of claim 20, wherein the said non-enhanced image and the said image produced after the contrast enhancement was infused are used to quantitatively measure changes of signal intensities in tissues caused by the contrast enhancement agent.
 22. A memory medium storing machine instructions for carrying out steps (b)-(d) of claim
 14. 23. A pulse sequence for quadruple inversion-recovery magnetic resonance imaging of a site, comprising: (a) a first double inversion pulse pair comprising in rapid succession and in either order a slice-selective inversion radiofrequency (RF) pulse and a non-selective inversion RF pulse; (b) a first time delay period; (c) a second double inversion pulse pair comprising in rapid succession and in either order a slice-selective inversion RF pulse and a non-selective inversion RF pulse; (d) a second delay period; and (e) an acquisition sequence comprising RF pulses and magnetic field gradient pulses for spatial encoding and readout of a signal for producing an image of the site.
 24. The pulse sequence of claim 23, wherein the slice-selective inversion RF pulse and the non-selective inversion RF pulse produce at least a partial if not a complete inversion of longitudinal magnetization.
 25. The pulse sequence of claim 23, wherein the first and the second delay periods are not equal.
 26. The pulse sequence of claim 23, wherein the first and the second delay periods are selected to ensure that a magnetization of blood is substantially zero when the acquisition sequence occurs.
 27. The pulse sequence of claim 23, further comprising more than two double inversion pulse pairs.
 28. The pulse sequence of claim 23, wherein said non-selective inversion pulses are rectangular pulses.
 29. The pulse sequence of claim 23, wherein said non-selective inversion pulses are adiabatic phase-modulated pulses.
 30. The pulse sequence of claim 23, wherein said slice-selective inversion pulses are adiabatic phase-modulated pulses.
 31. The pulse sequence of claim 23, wherein the said slice-selective pulses are targeted at one of a single slice, a group of slices, and a three-dimensional slab.
 32. The pulse sequence of claim 23, wherein the said acquisition sequence is used for one of two-dimensional and three-dimensional data acquisition.
 33. A system for carrying out quadruple inversion-recovery data acquisition, comprising: (a) an MRI apparatus adapted for producing an image of a site; and (b) a computer coupled to the MRI apparatus to control it, said computer including: (i) a memory in which machine instructions are stored; and (ii) a processor coupled to the memory, said processor executing the machine instructions to control the MRI apparatus to carry out a sequence of operations, including: (1) generating and applying in rapid succession a pair of radiofrequency (RF) inversion pulses, one RF inversion pulse being selective and the other being non-selective; (2) waiting for a first predefined delay period; (3) again generating and applying in rapid succession the pair of RF inversion pulses, one of which is selective and the other of which is non-selective; (4) waiting for a second predefined delay period; and (5) acquiring a magnetic resonance signal for use in producing an image of the site.
 34. The system of claim 33, wherein the machine instructions further cause the processor to determine the first and the second predefined delay periods so as to minimize a deviation of blood magnetization from zero, over a defined range of T₁ of blood at the site.
 35. The system of claim 33, wherein the machine instructions further cause the processor to determine the first and the second predefined delay periods so that blood magnetization at the site is substantially zero when the magnetic resonance signal is acquired.
 36. The system of claim 33, wherein the machine instructions enable an operator to input the first and the second predefined delay periods.
 37. The system of claim 33, wherein the machine instructions cause the processor to control the MRI apparatus so that the sequence is periodically repeated with a repetition time that is defined by an operator.
 38. The system of claim 33, wherein the machine instructions cause the processor to enable an operator to selectively run the sequence in synchronization with a cardiac rhythm of a subject that is being imaged, with a period equal to at least one systolic period. 